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ABSTRACT 

The IceCube Collaboration initially reported the detection of 37 extraterrestrial neutrinos in the TeV - PeV 
energy range. The reconstructed neutrino events were obtained during three consecutive years of data taking, 
from 2010 to 2013. Although these events have been discussed to have an extragalactic origin, they have not 
been correlated to any known source. Recently, the IceCube Collaboration reported a neutrino-induced muon 
event with energy of 2.6 =b 0.3 PeV which corresponds to the highest event ever detected. Neither the recon¬ 
structed direction of this event (J2000.0), detected on June 11 2014 at R.A.=110°.34, Dec.=ll°.48 matches 
with any familiar source. Long gamma-ray bursts (lGRBs) are usually associated with the core collapse of 
massive stars leading relativistic-collimated jets inside stars with high-energy neutrino production. These neu¬ 
trinos have been linked to the 37 events previously detected by IceCube experiment. In this work, we explore 
the conditions and values of parameters so that the highest neutrino recently detected could be generated by 
proton-photon and proton-hadron interactions at internal shocks inside 1GRB progenitor star and then detected 
in IceCube experiment. Considering that internal shocks take place in a relativistic collimated jet, whose (half) 
opening angle is 6q ^0.1, we found that lGRBs with total luminosity L < 10 48 erg/s and internal shocks 
on the surface of progenitors such as Wolf-Rayet (WR) and blue super giant (BSG) stars favor this multi-PeV 
neutrino production, although this neutrino could be associated to L ^ IO 50,5 10 50 ) erg/s provided that the 

internal shocks occur at ~ 10 9 IO 10,2 ) cm for a WR (BSG). 

Subject headings: Gamma-ray burst: general - stars: interiors - stars: jets - neutrinos 


1. INTRODUCTION 

The IceCube Collaboration has reported evidence for extragalactic neutrinos. They announced, firstly, two shower-like events 
of PeV neutrinos (Aartsen et al. 2013), after 24 events (19 shower-like, and the remainder track-like; IceCube Collaboration 
(2013)) and finally 11 events were added to these detections (Aartsen et al. 2014). In the aggregate, 37 events (9 track-like and 
28 shower-like) were detected after analyzing three year data (2010 to 2013; Aartsen et al. (2014)). The excess measured over 
the atmospheric background had a statistical significance 5.7cr. Arrival directions of these neutrino events seem to be isotropic 
and nonneutrino track-like has been associated with the location of known sources yet. Recently, a multi-PeV neutrino-induced 
muon event detected in 2014 June 11 was reported by the IceCube Collaboration (Schoenen & Raedel 2015). The neutrino with 
an energy 2.6 =b 0.3 PeV is the highest-energy event ever detected. This event had a probability of less than 0.01% of having an 
atmospheric origin, hence putting together in a skymap this event with the known sources as shown in Figure 1 , the reconstructed 
direction of this event (J2000.0: R.A.=110°.34, Dec.=ll°.48) does not match with any of the known sources. 

Long gamma-ray bursts (lGRBs) have been generally linked to core collapse of massive stars leading to relativistic-collimated 
jets and supernovae (CCSNe) of type Ib,c and II. Type II and lb are widely believed to have a radius of R* « 3 x 10 12 cm, 
usually associated with blue super giant (BSG) stars, and type Ic supemovae are thought to be Wolf-Rayet (WR) stars with radius 
R* ~ 10 11 cm. Depending on the luminosities and durations, successful lGRBs have revealed a variety of GRB populations: 
low-luminosity (11), ultra-long (ul) and high-luminosity (hi) GRBs (Meszaros & Waxman 2001; Liang et al. 2007; Gendre et al. 
2013). While llGRBs and ulGRBs have a typical duration of (~ 10 3 s), hlGRBs have a duration of tens to hundreds of seconds. 
Another relevant population associated with CCSNe, although unobservables in photons, are failed GRBs which could be much 
more frequent than successful ones, limited only by the ratio of type Ib/c and type II SNe to GRBs rates. This population has 
been identified by having high-luminosities, mildly relativistic jets and durations from some to ten seconds (Huang et al. 2002; 
Meszaros & Waxman 2001; Soderberg & et al. 2006, 2010). 

Many astrophysical environments have been explored to explain these high-energy (HE) neutrinos. For instance, GRB (Murase 
& Ioka 2013; Razzaque 2013; Fraija 2014; Liu & Wang 2013; Tamborra & Ando 2015a; Petropoulou et al. 2014; Tamborra 
& Ando 2015b), active galactic nuclei (AGN; Dermer et al. (2014); Stecker (2013); Petropoulou et al. (2015); Padovani et al. 
(2015)), dark matter (DM; Feldstein et al. (2013); Esmaili & Dario Serpico (2013)) and cosmic neutrino background connections 
(Ioka & Murase 2014; Ng & Beacom 2014). 

GRBs have been widely pointed out as emitter sources of HE neutrinos in connection to cosmic rays accelerated up to ultra 
high energies (Milgrom & Usov 1995; Waxman 1995b, a; Waxman & Bahcall 1999; Waxman 2000b, a; Vietri 1997; Rachen & 
Meszaros 1998; Bottcher & Dermer 1998). 

Recently, Murase & Ioka (2013) studied HE neutrino production in collimated jets inside progenitors of lGRBs, considering both 
collimation and internal shocks. Authors showed that whereas classical GRBs may be too powerful to generate HE neutrinos 
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inside stars, low-power (lp) GRBs (including llGRBs) were candidates for producing the HE neutrinos reported by the IceCube 
collaboration. 

One major progress in recent years has been the realization that photon emission from the photosphere, namely optically thick 
regions may be relevant to GRBs. Clearly, the photosphere marks the point below which the plasma is optically thick and above 
which it is optically thin. Several theoretical works in recent years dealt with the properties of shock waves that propagate in 
optically thick regions - the so called ’’radiation mediated shocks” (Budnik et al. 2010; Katz et al. 2010; Bromberg et al. 2011a; 
Levinson 2012). These authors have argued, that the properties of these radiation mediated shocks in the depth of the star are 
such that do not allow them to accelerate protons up to HEs. 

In this work, we explore the conditions under which the highest neutrino recently detected at 2.6 PeV could be created by hadronic 
interactions in internal shocks inside a GRB progenitor star. In section II, we show the dynamics of jet, internal shocks, particle 
acceleration process and the neutrino production by proton-photon (p- 7 ) and proton-hadron (p-h) interactions. In section III, we 
explain the neutrino flux expected on IceCube detector from p -7 and p-h interactions. In section IV we discuss our results and 
in section V we give a brief conclusions. We hereafter use primes (unprimes) to define the quantities in a comoving (observer) 
frame, Q x = Q/ ICE in c.g.s. units and k = h= c=l in natural units. 


2. INTERNAL SHOCKS AND NEUTRINO PRODUCTION 


In GRB context, one of the most prosperous theory to explain the prompt emission and the afterglow is the fireball model 
(Zhang & Meszaros 2004; Meszaros 2006). In this model, the prompt emission is explained through internal shell collisions, 
where faster shells (Tf) catch slower shells ( T s ). These internal shocks take place at a distance of 77 = 2T 2 t u , with t v given 
by the variability time scale of the central object and T cx yjY / T s the bulk Lorentz factor of the propagating shock. We are 
interested in those internal shocks that occur inside the star (r, < R* with R* the radius of the progenitor’s stellar surface) 
and the hydrodynamic jet is collimated by the cocoon pressure via collimation shocks (Bromberg et al. 2011b; Mizuta & Ioka 
2013). It happens when the jet luminosity is low and/or density is high for an initial opening angle. The total energy density 
U = 1/(8 7 r m v ) T- 4 Lt~ 2 in internal shocks is equipartitioned to amplify the magnetic field 


B' = e 1 / 2 r ~ 2 L 1/2 t 


-1 


( 1 ) 


with €b = Ub/U = (B 2 / 8 tt)U and accelerate particles ei = Ui/U , where i is the ith kind of particle with the condition 
cb + 6 i = 1. Electrons are accelerated in internal shocks and then are cooled down rapidly by synchrotron radiation in the 
presence of the magnetic field. The synchrotron photons in this length scale produce an opacity to Thomson scattering given by 
T ' th = T ~ 3 Lt~ x where gt is the Thompson cross section and m p is the proton mass. Photons can thermalize in a black 

body temperature with peak energy given by 


T' ~ — e ^/ 4 Li / 4 r -i +-i/2_ 


( 2 ) 


The proton distribution is cooled down in internal shocks via electromagnetic (synchrotron radiation and inverse Compton 
(IC) scattering) and hadronic p -7 and p-h interactions channels. P -7 and p-h interactions take place when accelerated pro¬ 
tons interact with thermal keV photons and proton density at the shock. In both interactions HE charged pions and kaons 
are produced; p + 7//1 -A X + it ± /K ± , and subsequently HE neutrinos 7r + -a /i + + -A e + + v e + and 

7r -A /i + -A e + E e + 


2 . 1 . Particle acceleration process 

Before the jet breaks out the star, it keeps the bulk Lorentz factor low by converging cylindrically via collimation shocks under 
the cocoon pressure. The characteristic radius to which the shocked jet becomes cylindrical is (Bromberg et al. 2011b) 
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where p a is the ambient density, 60 is initial opening angle, £ c = 3777 , Vc — 3 _ a§ _ 7 _ as ■ 

absolute jet luminosity defined through the total luminosity L. In the meantime, the position of the non-relativistic head is 
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with K given by ( 7 r ) 17 ^ 7 for r a < r < and (t ^) 17 ^ 7 / r _^*l 5 for r > r^. Here po = 3.4 x 10 5 gem 3 , r a = 10 10,8 


(■ r b -R *) 5 


cm, 


r\) = lO 1 ^ cm, M * is the progenitor mass and as = 3/2 - 3, being 3/2 and 3 for convective and radiative envelopes, respectively. 
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From eqs. (3), (4) and (5), the cylindrical and head radii for a WR are 


and for a BSG are 
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Propagating shocks in optically thick regions are radiation mediated. The properties of these radiation mediated shocks in the 
depth of the star are such that do not allow to accelerate protons up to HEs. Upstream proton flow is decelerated by photons 
created in the downstream and diffusing into the upstream region (Budnik et al. 2010; Katz et al. 2010; Bromberg et al. 2011a; 
Levinson 2012). It happens when the deceleration scale Idee ~ (n p cr T y±)~ 1 of the shock width is smaller than the comoving 
size (l v ) of the upstream flow, as indicated by Murase & Ioka (2013). Protons, then, are expected when shock propagates in an 
optically thin region (radiation unmediated shock; l p <C Idee )• Hence, accelerated protons are expected around the photosphere 
r ~ (1 — 10) (Murase 2008), as assumed in the dissipative photosphere scenario (Rees & Meszaros 2005; Murase et al. 2006). 
Following Murase & Ioka (2013), a reasonably necessary condition for proton acceleration is 

r = n p arr < min[r 2 eZ ,0.1C _1 r 3 eZ ] , (8) 

where C = 1 + 2 In T 2 rel , T rei « (I7/r + Y/Y f)/2 is the relative Lorentz factor between the faster and merged shells and n p is 
the number proton density. 

In this case, neutrino energy is limited by the maximum energy at accelerated protons which is obtained equaling the time scales 
of acceleration t' acc = 27r ^ e,p E' p e^ 1 ^ 2 L -1 / 2 T 2 t v and cooling synchrotron t' p syn = L~ x T 4 £ 2 , supposing that 

the synchrotron process is the longest timescale (e.g. see Fraija (2014)). Hence, it can be written as 
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where e is the electric charge, m e is the electron mass, £ is a factor of equality and B' is the critical magnetic field for protons. 


2.2. P -7 interactions 

The number density of thermalized photons t]' ~ 2< ^P P can be explicitly written as 


n' = — ' ) e 3 J A L 3/i r~ 3 t~ 3/2 . (10) 

The optical depth characteristic for this process is r' = 6-5 T -2 ^ 1 ^ 2 , where ((x) is the zeta function and a P1 

is the cross section of p -7 interactions. From the resonance p + 7 —» A + in the rest frame and taking into account that each 
neutrino carries 5% of the proton energies (E v = 1/20 E p ), then the threshold neutrino energy can be written as 


where m/± is the resonance mass. 
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2.3. P-h interactions 


Protons co-accelerated in the internal shocks have a comoving number density (Razzaque & Smirnov 2010; Razzaque et al. 
2005; Ando & Beacom 2005; MacLachlan et al. 2013) 


n p = 1/(8 7 rm p ) T 4 Lt u 2 , 


( 12 ) 


with energy density n' p m p . The optical depth for p-h interactions is r pp ~ 47777 LV 3 t v 1 where a pp = 34.3 + 1.88 S + 

0.25 S 2 mb is the cross section for these interactions with S = log (E p /TeV) (Kelner et al. 2006). Pions ( 7 r ± ) and kaons (Ff ± ) 
can be produced as products of p-h interactions (p + h -» X + 7 r^/A^). Considering the meson lifetime and p-h interaction time 
scale 


4 ' 107rm p r _i-p4 + 2 

7 had ~ “ ^ 1 l v j 


J PP 


(13) 


mesons, then, in the presence of the magnetic field radiate via synchrotron radiation and also could interact with the primary 
protons before they decay. 
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2.3.1. Pion component 


O'/i 11 L Hj 

Equaling the cooling synchrotron radiation t' n+ syn = ^ ^ e^ 1 L~ x T 2 E'~^ and the lifetime of pions t' n+ dec = 


e' _ 


7 r+,dec m _|_ 

with the hadronic time scale (eq. 13), we get neutrino break energies due to the synchrotron radiation and the lifetime of pions 
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respectively. Here ra^ is the charged pion mass and r 7r + is the lifetime of the charged pion. 


2.3.2. Kaon component 


Comparing for kaons, the cooling synchrotron radiation t' k+ syn = L 1 T 2 t 2 E + 1 and the lifetime t' k+ dec = §^r k + 

with the hadronic time scale (eq. 13), we obtain neutrino break energies due to the synchrotron radiation and the lifetime of 
kaons 
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respectively. Here m k + is the charged pion mass and r k + is the lifetime of the charged kaon. 


3. NEUTRINO EXPECTATION 

The expected number of reconstructed neutrino events in the IceCube experiment can be computed as 

N ev = TN A f a uN M e ff —A dE v , (18) 

J E* h 

whereT cx 4 years is the observation time, NU=6.022x 10 23 g _1 is the Avogadro number, cr v n = 6.78x 10 _ 35 cm 2 (£ T zy /Tey ) 0,363 
is the neutrino-nucleon cross section (Gandhi et al. 1998) and M e ff is the effective target mass of the IceCube experiment (Ice- 
Cube Collaboration 2013). The neutrino spectrum dN v /dE v is computed from p -7 and p-h interactions as follows. 


3.1. P -7 interactions 

The spectrum of muon neutrino generated by p -7 interactions is (Kelner & Aharonian 2008) 
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where V = f ?rr 3 , d z is the luminosity distance from the source and the proton distribution is given by a power law 

dN P ~ £ fi-OL 

dE~ Aph * 5 


(19) 


( 20 ) 


with Ap ~ n p /GeV. The distribution of target photons dN e /de corresponds to the blackbody radiation with temperature and 
photon density given by eqs. (2) and (10), respectively, and the function F^(rj, x) is given in appendix A. 


3.2. P-h interactions 

The spectrum of muon neutrino generated by p-h interactions is (Kelner et al. 2006) 
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where the proton distribution dN p /dE p is given in eq. (20). The function F v = F ( i> + F ( 2 ) comes from the contribution of 
decay of muon F ( 1 ) and from the direct decays of pion F ( 2 ). Functions F ( 1 ) and F ( 2 ) are explicitly written and explained in 

V p Vp V p ISp 

appendix B. 
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4. RESULTS 

In the CCSNe-GRB connection framework, we have considered WR and bigger progenitors like BSG stars with formation of 
jets leading to internal shocks inside of them. In these shocks, energy is equipartitioned between magnetic fields and acceleration 
of particles. Electrons and protons are expected to be accelerated in these shocks and after, to be cooled down by synchrotron 
radiation, inverse Compton and hadronic processes (p -7 and p-h interactions). Both thermalized photons produced by electron 
synchrotron radiation and hadrons at the shocks serve as targets in the neutrino production through hadronic interactions. In this 
model, we have shown the production channels of HE neutrinos as a function of bulk Lorentz factor, variability time scale, total 
luminosity and microphysical parameters (e# and e e ). A detailed description such as time scales for all processes (electromagnetic 
and hadronic) is given in Fraija ( 2014 ). Additionally, we have computed the number of events expected in the IceCube experiment. 
In Figure 2 is shown the regions for which internal shocks occur inside a WR (red region) and BSG (blue region) star. By 
considering the typical values of variability 10 _3 < t v < 1 s, we have found that the bulk Lorentz factor lies in the range 
1 < T < 40 for a WR and 1 < T < 220 for a BSG. 

From eqs. ( 6 ) and (7), and using the typical values of #0 and L (Mizuta & Ioka 2013 ), the cylindrical and head radii for a WR 
and BSG can be written as 


1.01 x 10 8 cm t{- 08 Lgi 81 0g;I?! i£; 11 for WR 

3.38 x 10 9 cm t% /5 L % 1 8 ° 0^ for BSG , 


( 22 ) 


and 


9.67 x 10 10 cm f}' 03 L% i 35 ^ 9 ® 9 M" 9 ^ forWR 

2.29 x 10 12 cm t\ /6 L % 1 ° 7-i for BSG . 


(23) 


Here we have used a$= 2.1. Taking into account the condition r c < rj and eq. ( 8 ), we plot contour lines of total luminosity and 
internal shock radius as a function of bulk Lorentz factor for which the necessary condition for proton acceleration is established 
in a WR (left panel) and BSG (right panel), as shown in Figure 3. From the left panel (WR), one can see that the total luminosity 
lies in the range 10 46 ' 5 < L < 10 50 ' 5 erg/s for an internal shock radius in the range 10 8 < Tj < 10 11 cm and from the right 
panel (BSG), the total luminosity lies in range 10 46 < L < 10 49,8 erg/s for a radius 10 10 < rj < lO 12,3 cm. When we consider 
the maximum values of bulk Lorentz factors T ^ 40 and T ^ 220 for a WR and SGB (see fig. 2), respectively, from eqs. ( 8 ) and 
( 12 ), the optical depths become 


f 2.11 £(50.5-48)7/ 


(24) 


Comparing eqs. (22) and (24), it is trivial to demonstrate that accelerated protons could be expected when r c < rj and total 
luminosity in the range 10 48 < L < 10 50,5 erg/s (10 48 < L < 10 5 ° erg/s) for a WR (BSG). For instance, protons accelerated on 
the surface of either a WR or BSG may be expected under the condition that the jet would have a luminosity Lj < 10 48 erg/s for 
both a WR and BSG. Otherwise, under the condition r 3 < r c protons cannot be expected inside the progenitors unless L > 10 51 
erg/s or/and t v < 10 -3 s and T > 40 (T > 220) for a WR (BSG). Comparing the acceleration and the synchrotron time scales, 
then the maximum energy achieved by accelerated protons in internal shocks with r c < rj and r ~ 1 is given by eq. (9). 

Figure 4 shows neutrino energy (above 30 TeV) generated by p -7 interactions at internal shocks inside of a WR (left; n 7 = 
10 16 cm -3 ) and BSG (right; n 7 = 10 13 cm -3 ) for a luminosity in the range 10 46 < L < 10 52 erg/s and four values of energy 
fraction going to accelerate electrons (e e =0.5, 5 x 10 -2 , 5 x 10 -3 and 5 x 10 -4 ). Left panels exhibit that the maximum 
neutrino energy for L = 10 52 erg/s and e e = 0.5(5 x 10 -4 ) is E v = 41.7(9.8) x 10 15 eV whereas the neutrino energy is 
E v — 5.3(30.1) x 10 16 for L = 10 46 erg/s and the same values considered of e e . As shown in the right panels, the maximum 
neutrino energy is E v = 1.4(8.1) x 10 14 eV for L = 10 52 erg/s and e e = 0.5(5 x 10" 4 ), and E v = 4.2(25.3) x 10 15 eV for 
L = 10 46 erg/s. In panels can be observed a dashed line with E u = 2.6 PeV. This value of neutrino event can be obtained for 
L < 7.3 x 10 46 (7.1 x 10 49 ) erg/s and e e < 0.5 (5 x 10 -4 ) when considered a WR, whereas it is created for L < 10 51 erg/s 
and e e < 0.5 in a BSG. Finally, based on the condition required for the proton ray acceleration (eq. 24), we have highlighted the 
region where parameters are restricted, as indicated in the plot labels. 

Figure 5 shows neutrino energy (above 30 TeV) generated as charged kaon and pion decays in p-h interactions for 10' - 4 < 
e B < lir 1 (left) and 10 46 < L < 10 52 erg/s (right). Left panel shows that the maximum neutrino energy is E v < 3 x 10 13 
(2.1 x 10 16 ) eV for e# = 10 -1 and E v — 7.1 x 10 14 (> 10 19 ) eV for cb = 10 -4 . Right panel exhibits that the maximum 
neutrino energy is E v = 1.1. x 10 17 (4.5 x 10 21 ) eV for L = 10 45 erg/s and E v < 3 x 10 13 (5.8 x 10 14 ) eV for L = 10 52 erg/s 
when neutrino is coming from pion (kaon) decay product inside WR (BSG). Putting a dashed line to mark the neutrino energy 
E v - 2.6 PeV, we can observe that for cb < 10 -3 (5.1 x 10 -3 ) neutrinos at this energy could be generated from decay product of 
pions (kaons) in a WR whereas for 6.1 x 10 -3 < e# < 0.1 it could be created in a BSG. The right panel shows that depending 
on total luminosity, neutrinos coming from a WR and/or BSG could be created as pion and kaon decay products. In this panel, 
we show the region where parameters are restricted as showed in the plot label in accordance with the condition required for the 
proton acceleration (eq. 24). 

Requiring the effective target mass of the IceCube experiment (IceCube Collaboration 2013), we use the method of Chi-square 
X 2 minimization as implemented in the ROOT software package (Brun & Rademakers 1997) to fit the effective target mass as a 
function of neutrino energy. Then, the best function that describes the effective target mass is 


f f 6 (E v ) , 1 TeV < E v < 950TeV 

\ 5.19 x 10 “ 3 (^ 7 ) + 3.86 x 10 2 , 950TeV < E v < 10 4 TeV, 
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It is worth noting that the muon neutrino of the effective target mass was used. Figure 6 shows the effective target mass of the 
IceCube experiment as a function of neutrino energy. 

With the p -7 spectrum (eq. 19) and using the method of Chi-square x 2 minimization (Brun & Rademakers 1997), we have 
obtained the functions (eqs. A3, A4 and A5) that describe the parameters a\, a 2 and as, respectively, as shown in Figure 7 and 
Appendix A. 

Taking into account the target mass of the IceCube experiment (eq. 25), the p -7 (eq. 19) and p-hadron (eq. 21) spectra, from 
eq. (18) we have computed the number of neutrinos expected in the IceCube detector. Figure 8 shows the number of neutrinos 
as a function of energy expected in the IceCube experiment from the p -7 (above panels) and p-h (below panels) interactions 
for a source located at 10 Mpc. From p -7 interactions we have used a target photon density n 7 ^ 10 13 cm -3 whereas from 
p-h interactions a proton density target (as indicated in the label plots) was required. We compute the upper limit on the proton 
density in order to obtain one neutrino event at 2.6 PeV for typical values of the spectrum index of proton distribution (2.1 < a < 
2.3). 

5. CONCLUSIONS 


In order to explain the multi-PeV neutrino, we have evoked p -7 and p-h interaction model at internal shocks inside the 1GRB 
progenitor stars. We have explored the set of parameter values: T, t u ,L, a, cb and e e to explain the neutrino event at 2.6 PeV. 
We have computed the characteristic radius and the optical depth to which the shocked jet becomes cylindrical and the condi¬ 
tion for proton acceleration inside the star is satisfied. We have found that protons can hardly be accelerated at internal shocks 
formed at the distance less than the cylindrical radius 77 < r c unless L > 10 51 erg/s or t v < 10 _3 s and T > 40 (T > 220) 
for a WR (BSG). Otherwise, protons at internal shocks under the condition r c < 77 can be accelerated, thus interacting with 
hadrons and photons at this place. We have found that this neutrino could be created by p -7 and p-h interactions with a total 
luminosity and internal shock radius in the range 10 46,5 < L < IO 50,5 erg/s (10 46 < L < 10 49,8 erg/s) and 10 8 < 77* < 10 11 
cm (10 10 < 77 < IO 12,3 cm), respectively, for a WR (BSG) with T < 40 (T < 220). Additionally, the values of equipartition 
parameters lie in the ranges e e < 0.05 and 10 -4 < €b < 10 _1 which could be a smaller range depending on the progenitor 
considered. 

Recently, Murase & Ioka (2013) studied GeV - PeV neutrino production in collimated jets inside progenitors of GRBs, consid¬ 
ering both collimation and internal shocks. They showed that although internal shocks in the collimated jet are less favorable 
for the cosmic ray acceleration than the collimation shocks, due to the meson radiative cooling HE neutrinos around 1 PeV can 
be only expected from internal shocks. Otherwise, HE neutrinos from the collimation shocks are expected around 10 TeV due 
to the strong meson cooling. They also found that unless T > 10 3 , cosmic rays and HE neutrinos are expected for high-power 
jets inside a WR. Authors showed that llGRBs can be consistent with the astrophysical neutrino background reported by IceCube 
Collaboration. In this work, we have explored the values of total luminosity, internal shock radius, bulk Lorentz factor and the 
microphysical parameters so that CCSNe-GRB located at 10 Mpc (Ando & Beacom 2005) could create a 2.6-PeV neutrino event 
at internal shocks inside the star. 

In hlGRBs copious target photons are expected for photo-hadronic interactions, making these GRBs promising candidates for 
neutrino detection. Searches for HE neutrinos in spatial and temporal coincidence around this population has been performed, 
although no neutrinos have been observed. The null neutrino result reported by IceCube Collaboration in some GRBs (Abbasi 
et al. 2012) could be explained based on the results obtained and showed in figs. (3), (4) and (5). We show that a multi-PeV 
neutrino event could be associated to lGRBs with L ~ IO 50,5 (^ 10 5 °) erg/s provided that the internal shocks occur at ~ 10 9 
(~ IO 10,2 ) cm for a WR (BSG). Otherwise, this PeV neutrino is not expected from hlGRBs with T <C 40 (T <C 220) and/or 
internal shocks on the surface of WR (BSG). 

By considering a range of GeV - TeV neutrino energies with a flux ratio of N v cx ~ 2N Ue ~ 2 N De , Fraija (2015) studied 
the effect of thermal and magnetized plasma generated in internal shocks and the envelop of star on the neutrino oscillations. The 
author showed that the resonance lengths lie in the range l res ~ ( 10 12 to 10 14 ) cm, hence depending on the progenitor associated 
(e.g. WR or SBG) and oscillation parameters, neutrinos would leave the internal shock region and the star in a different flavor 
to the original. Fraija (2015) found that although neutrinos from a few GeV to TeV can oscillate resonantly from one flavor 
to another, neutrino > PeV can hardly oscillate. Finally, he showed that a non-significant deviation of the standard flavor ratio 
(1:1:1) could be expected on Earth. 

The number of sources with internal shocks inside the stars may be much larger than those exhibiting one. Within 10 Mpc, the 
rate of core-collapse supernovae is ~1 - 3 yr -1 , with a large contribution of galaxies around 3-4 Mpc. At larger distances, the 
expected number of neutrino events in the IceCube detector is still several, and the supernova rate could be > 10 yr -1 at 20 Mpc 
(Ando & Beacom 2005). Interference effects in the detector by atmospheric neutrino oscillation are very small (less than 10 %) 
due to the short path traveled by neutrinos in comparison with cosmological distances (Mena et al. 2007). 


We thank C. Dermer, A. M. Sodelberg, B. Zhang, I. Taboada, K. Murase, W. H. Lee and F. de Colle for useful discussions, 
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Fig. 1.— Sky-map with the location of BL Lac sources, FSRQs, Supernova Remnants, Radio Galaxies and the multi-PeV neutrino event detected by the 
IceCube detector (blue cross). 

Anita Miicke and Antonio Galvan for helping us to use the SOPHIA program and the TOPCAT team for the useful sky-map tools. 
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APPENDIX 
A. FUNCTION F t 


P 7 


The function F P1 ( 77 , x) given in eq. (19) is (Kelner & Aharonian 2008) 

F^(r],x)=o„ exp j -a 2 log (jQ J log 


where ib = 


ip = 2.5 + 1.4 log with 77 =s 


_ J-jp 


and 770 = 2 ^ + ^ ~ 0.313. The variables y and x are 

m P m p 

x — x\ . E v 

— - and x — — , 

E v 


(Al) 


X2 — Xi 


and 


(A2) 
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Fig. 2.— Set of values of variability and bulk Lorentz factor for which internal shocks take places inside a WR (red region) and BSG (blue region). 




r-j (cm) 


r-j (cm) 


Fig. 3.— Necessary condition for proton acceleration (eq. 8) when we consider a WR (left) and BSG (right). 

respectively. The appropriate parameters of the distribution for muon neutrino are: x\ = 0.427.x_ and 

( 0.427 x + , 77/770 <2.14 

X2 = < [0.427 + 0 . 00729 ( 77/770 - 2.14)] x+ , 2.14 < 77/770 < 10 

l £+ , 77/770 < 10. 

The values of x± correspond to the maximum and minimum pion energies, x± = 2 (i+ri) fa + r2 =*= ^(77 - r 2 - 2 r )(77 - r 2 + 2 r)\ 
with r — m^/rrip « 0.146. 

In accordance with the values a±, a 2 and a 3 given in Kelner & Aharonian (2008)) and using the method of Chi-square x 2 
minimization (Brun & Rademakers 1997), the best fit of the numerical values of parameters a±, <22 and a 3 are 

01 = 1.63 x 10“ 18 (—) + 1.49 x 10 -16 (-)- 1.65 x 10“ 16 , (A3) 

\Vo) \VoJ 


f 6 (n) ’ 1< r//r/o < 9 

1.94 X 10- 5 (^) 2 - 2.87 x 10- 3 (i) + 0.14 77/770 >9. 












(eV) ^ E v (eV) E v (eV) 
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L (erg/s) L (erg/s) 



.4.— Neutrino energy as a function of luminosity generated by p -7 interactions at internal shocks inside of a WR (left) and BSGs (right). 
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Fig. 5.— Neutrino energy generated by p-h interactions at internal shocks, synchrotron (left) and lifetime (right). 
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Neutrino Energy (TeV) 

Fig. 6 .— Effective target mass of the IceCube experiment as a function of neutrino energy (green line; IceCube Collaboration (2013)). Line in blue color 
represents the function used to fit the effective target mass. 





Fig. 7.— Fit of numerical values of parameters a±, a 2 and <23 characterizing the neutrino spectrum (Kelner & Aharonian 2008). 
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r_j = 10 l:L cm 


-> in 12 
rj=3 x 10 cm 




Fig. 8 .— Number of events as a function of neutrino energy expected in the IceCube experiment from the p -7 (above panels) and p-h (below panels) interactions. 


with 


and 


with 



&3 


’ 1 < 77/770 <4 

- 1.74 x 10- 4 (^) 2 + 2.25 x IQ ” 2 (^) + 2.10 77 / 77 0 > 4. 


U 



ss —0.97 



+ 9.81 



- 3.52 x 10 


+53.32 



- 27.81 



(A4) 


(A5) 


respectively. 
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B. FUNCTION Ff p 

The function F pp is split in two parts (Kelner et al. 2006): one coming from the pion decay 7 r uis u (F p F) and the other 

comes from decay of muon /i —( F pp 2) ). 


Bl. F pp r) from pion decay 

The spectrum of muonic neutrino generated through de direct decay of pion can be written as 


Fm = 61 


log z f 1 - z b2 
z \1 + b 3 z b *(l - z b *) 

1 4 b 2 z b2 Ab 3 b 2 z b2 (l-2z b2 ) 


_log z l-z b 2 1 + b 3 z b2 (l - z b *) \ ’ 

where the variables z and x are z = and x = ff, respectively. The parameters b\, b 2 and 63 are 

61 = 1.75 +0.204 S + 0.010 S 2 , 

1 


b 2 = 


and 

with S = log (E p /TeV). 


1.67 + 0.111 S +0.0038 S 2 ’ 
b 3 = 1.07 - 0.086 S + 0.002 S 2 , 


(Bl) 


(B2) 

(B3) 

(B4) 


B2. F P (! 2) from muon decay 

The spectrum of muonic neutrino generated through de direct decay of muon can be written as 

[1 + c 3 (loga ;) 2 ] 3 


where the parameters c \, C 2 and C 3 are 


and 


F(2) = C\~ 


c 1 = 


1 


-(-log a ;) 5 


69.5 + 2.65 S +0.3 S 2 
1 


C2 = 


(0.201 + 0.0062 S + 0.00042 S 2 ) 1 ^ ’ 
0.279 +0.141 S +0.0172 S 2 


C3 = 


0.3 + (2.3 + S) 2 


(B5) 

(B 6 ) 

(B7) 

(B 8 ) 














